arXiv:hep-th/9701098vl 20 Jan 1997 


SOURCE DRIVEN SOLUTIONS 
OF QUANTUM FIELD THEORIES^ 

G.S. GURALNIKg 


Center for Theoretical Physics 
Laboratory for Nuclear Science 
and Department of Physics 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 


1 Introduction 

This talk describes a program to develop a new numerical approach to the solution of 
quantum field theories |j|] which we call the Source Galerkin method. The ideas involved 
have little directly in common with the Monte Garlo techniques that have demonstrated 
the power of numerical approaches to obtain results from quantum field theories. Source 
Galerkin has strengths which should allow its application to problems, such as scattering 
predictions, that require too much computing power to be examined by other methods at 
this time. Unlike Monte Carlo evaluations, Source Galerkin is equally suited for studying 
fermionic or bosonic problems. Further, at least in lower dimensions, a continuum formula¬ 
tion of Source Galerkin is straightforward. We have discovered that, while Source Galerkin 
implementations appear to be relatively economical in their use of computer cycles, they 
require complicated algebraic set up procedures for each problem examined. These take 
a good deal of care and large amounts of human time. While we have made considerable 
progress in understanding the issues in defining and applying these ideas, we have not yet 
been able to complete enough calculations to evaluate the chances that they might be able 
to supersede Monte Carlo techniques in any of their traditional applications. We have, 
however, already obtained a striking result arising from our ideas. This result consists of a 
deep new insight into the solutions of QFT and phase structure |§, H], |^. As one aspect of 
this, we have shown that the conventional definition of the path integral does not include 
all solutions of canonical quantum field theory because a very restricted set of boundary 
conditions is implicitly assumed. The full solution set can be obtained by appropriate con¬ 
tinuations of the integration regions into the complex plane. This insight was a result of the 
very detailed and specific understanding of a theory that is necessary in order to numerically 
perform evaluations using Source Galerkin methods. 
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To apply the Source Galerkin ideas, we first write a field theory action in the presence of 
external sources for the fields. Then, by enlightened guessing, an approximate set of solu¬ 
tions to the linear source functional differential equations for the vacuum-vacuum amplitude 
Z of this action is constructed. In principle, a terrible initial guess could be made. Through 
iteration this usually will eventually converge to a correct solution. For a QFT it is possible 
and desirable to utilize symmetries expected of the exact answer in the construction of the 
approximate solutions. This possibility is one of the great advantages of our method. In 
contrast it is difficult, or impossible, to insert knowledge obtained from analytic techniques 
into the construction of solutions using Monte Carlo techniques. 

A straightforward way to make a Galerkin guess is to pick some linear combination 
of functions from a subset of a complete set of functions of the source variables for the 
system being studied. The use of a complete set of functions has the advantage of a 
clear iteration path defined through increasing the number of functions used. Generally 
a guess will not satisfy the original equations. However, by picking the coefficients of the 
linear combination correctly, it is possible to pick an optimal approximation to the correct 
answer. The idea of Galerkin is to fix these coefficients by adjusting them so that the guess 
satisfies the differential equations as a weighted average over some region of source space 
1^] . The optimal definition of average is part of the problem of defining a rapidly convergent 
approximation scheme. Generally, the functions used to weight the average are chosen from 
the same set of complete functions as were used to build the guess. Gareful implementation 
of this scheme allows the reduction of the error between the approximate solution and the 
exact answer as the number of iterations is increased. 

A crude, but very illustrative, way to calculate a field theory is to expand .E in a Taylor 
series in the sources and then truncate this series at some maximal power i I 0- Such 
truncated series are dehnitely not solutions to non-trivial field theories. Nevertheless, it 
is reasonable to hope, for small sources, that the contributions from higher terms in the 
expansion rapidly decrease. We have checked and confirmed this on many problems where 
we know the answer or can check our results using Monte Carlo. The weight terms for the 
averaging procedure were usually chosen from the derivatives of the polynomial set used to 
approximate Z. The accuracy that is achieved with a few terms is amazing and holds over 
a wide range of parameters. 

Galerkin approximation is a weighted residual procedure. Such methods guarantee 
convergence in a mean. Our procedure is constructed so that the integral of the square of 
the local error vanishes in an appropriate limit. The more familiar variational Hamiltonian 
approximation techniques are in the same general class of approximations. However, Source 
Galerkin methods do not have a direct correspondence to such methods as they have been 
used in high energy physics. 

2 Example application of Galerkin Method 

Before we examine a field theory, we outline the application of Source Galerkin to a generic 
linear differential equation 


L/ = 0. 


(1) 


Here L is a differential operator defined in some domain D with boundary conditions defined 
on a surface dD. We begin by guessing a solution 

M 

f* = ipo{xi,... ,xn) + ’®^)- (2) 

i=i 

We assume the (pj are linearly independent and that they satisfy the appropriate boundary 
conditions but are otherwise unrestricted. 

Now make the definition: 

M 

R = Lr = + ( 3 ) 

i=i 

Next, pick M linearly independent test functions pi (xi, ... ,X]\f). We do not confine our con¬ 
siderations to the usual Galerkin choice of picking pi = (pi. Fix the coefficients ai,... ,aM 
by requiring 

(i?,p,) = 0, i = (4) 

We dehne this generalized dot product defined as follows; 

{R,Pi) = j dxi .. .dxNR{xi,... ,XN)pi{xi,... ,XN)d{xi,... ,xn)- (5) 

To make sure that this is sensible, d is picked so that it dampens integrals. Options that 
we have used in field theory problems include restricting d to be composed of 5 functions 
or choosing it to be of the form exp — XiAijXj . After choosing d and performing the 
above integrations, we obtain M linear equations in M unknowns. If these are independent, 
we can solve for all of the unknowns. If not, we must change our choice of test functions 
or increase their number. Under very general circumstances these approximations converge 
(weakly) to the exact answer as M ^ oo. 

3 Fermions 

An exciting feature of our application of Galerkin methods to field theory is that fermions 
can be included in a natural manner. Fermions present difficulty in Monte Carlo approaches 
for two reasons. The first is because of the “fermion determinant”. The problem occurs 
because after an approximate bosonic configuration is generated, the fermionic information 
contained in a determinant involving every point of the lattice system must be evaluated. 
This evaluation requires a large amount of computer time and it must be done for many 
bosonic configurations. There are, of course, many schemes to decrease the amount of 
computation this requires. However, the difficulty is intrinsic to how path integration forces 
us to deal with fermions. Source Galerkin avoids this difficulty because it allows making 
approximations for fermions without first solving a bosonic problem. I will outline how 
this is done in the following. The second problem is the well known “fermion doubling 
problem”. This difficulty has nothing to do with numerical methods, but is an artifact 



of lattice formulation. As mentioned, it is possible that we can avoid this as well. We 
have worked several fermionic problems in low dimensions on the continuum. This has 
been straightforward because the divergences we dealt with were trivial. We are currently 
working on a “proof” that Galerkin has enough flexibility to regularize all calculations 
without being forced to resort to a lattice. 

The key difficulty in developing a Source Galerkin procedure for fermions is that the 
naive inner product definition is not appropriate because of the Grassmann nature of the 
sources. Direct integration of Grassmann polynomials is defined through the usual integra¬ 
tion rules leads to a decimation of fermionic polynomials. This makes it difficult to generate 
a formula to analyze fermionic problems in an easy symmetric manner. To avoid this, we 
define a similar inner product to that developed for bosons. The definition we use is 

{A[fi,7]],B[fj,r]]) = J[dfi][dr]]exp[-fi7]]A[fj,r]]B[fi,7]]. (6) 


4 Numerical Solutions 


To test the source Galerkin ideas as applied to fermionic systems, we have studied simple 
quartic interactions on a one dimensional lattice. These are described by the action; 
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This action leads to the source functional differential equation: 
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If we restrict our study to 4 sites, this action can be solved exactly because of natural 
truncation of the power series expansion in the external sources at eighth order. Because 
we know the solutions, this model provides a non-trivial test case to examine the validity of 
our Galerkin methods. We perform this test by assuming approximate solutions truncated 
at fourth and then sixth order in the fermionic sources. The weights of the polynomials 
in the sources are set by the Galerkin procedures outline above. To see exactly how these 
answers converge, we set M = 1 and substitute a range of values for G and examine the 
behavior of the coefficients. We find that the convergence is rather good. This is shown in 
Tables 1 thru 3 where we display the values for several couplings for the three independent 
values of the two-point Green’s function. 

While the 4 site problem is relatively straightforward, even a five site lattice is beyond 
simple direct solution. However, the Galerkin results can be extended to a lattice of any size 
as long as the approximating polynomial is kept to lengths that are possible to handle in 



g 

4th Order 

6th Order 

Exact 

0 

0.75 

0.75 

0.75 

0.1 

0.7003 

0.7089 

0.7095 

0.5 

0.5553 

0.5840 

0.5926 

1.0 

0.4428 

0.4803 

0.5 

10 

0.0996 

0.1160 

0.1494 


Table 1; Convergence of fermionic Galerkin method for 


g 

4th Order 

6th Order 

Exact 

0 

0.25 

0.25 

0.25 

0.1 

0.2413 

0.2334 

0.2328 

0.5 

0.2130 

0.1928 

0.1852 

1.0 

0.1871 

0.1677 

0.15 

10 

0.0626 

0.1141 

0.0390 


Table 2; Convergence of fermionic Galerkin method for {ip(2)ip{l)) 


g 

4th Order 

6th Order 

Exact 

0 

0.25 

0.25 

0.25 

0.1 

0.2299 

0.2223 

0.2217 

0.5 

0.1726 

0.1543 

0.1481 

1.0 

0.1299 

0.1118 

0.1 

10 

0.0201 

0.0190 

0.0065 


Table 3; Convergence of fermionic Galerkin method for {ip(3)ip{l)) 























4th Order 

6th Order 

Mean Field 

0.5458 

0.5831 

0.6287 

0.2288 

0.2104 

0.2218 

0.0883 

0.0753 

0.0797 

0.0177 

0.0329 

0.0243 

0.0166 

0.0267 

0.0196 


Table 4: Interacting Propagator on 8 Sites with M = 1.0, g = 0.5. 

a computer. This is still a considerable restriction since the number of source polynomials 
eventually grows exponentially with the number of sites on the lattice. Nevertheless, if the 
higher order Green’s functions of the theory are evenly moderately well behaved, we would 
expect an even limited expansion combined with Galerkin averaging to give a good answer 
for the lower order Green’s functions. We have calculated the two point propagator for a 
eight site lattice using a Grassmann polynomial truncated first after fourth and then after 
sixth order. This is a non-trivial calculation since the fourth order polynomial by itself has 
has 124 coefficien ts. All coefficients were set by the Galerkin method. These results are 
displayed in Table 4. 

In Table 4, we compare our results with fourth order and sixth order Galerkin on an 
eight site lattice to those obtained by a mean field calculation. Mean field theory is a single 
pole approximation to the propagator. The Source Galerkin calculated propagator, on the 
other hand, should include contributions from all poles. Nevertheless, for this problem the 
answers are very similar. Calculations with the Source Galerkin method are very efficient, 
especially when compared to Monte Carlo. Usually, for lattice gauge theory, the fermions 
must be quenched to make calculations tractable. Here, determination of the interacting 
propagator for a system of dynamical fermions presents no special difficulty even in a pure 
fermionic theory. The bulk of a calculation involves a single matrix inversion for a given set 
of parameters. This is in contrast to Monte Carlo where many sweeps through the lattice 
are necessary to reduce statistical error. As our tables illustrate, the Source Galerkin 
calculations are very clean and show rapid convergence, when compared to the mean field 
results even at intermediate couplings and using only low order polynomials. 


5 Non-Linear Methods 

Our initial examinations of the Galerkin method concentrated on expansions that were 
simple for Z. However, because of the disconnected Green’s functions content of Z, it is 
often much better to deal with In Z. The difficulty with this is that the resulting equations 
are non-linear in the Galerkin parameters. For example, from the previous equations, it is 
straightforward to show that the approximation for a solution to quartic scalar field theory 
of the form 


.E* = e 


Jj-\- JlJm-^lmpqJpJq'^ 


(9) 




yields non-linear equations after being acted on by the source differential equations and then 
Galerkin averaged. These can be handled numerically. It is possible to develop iteration 
schemes based on the above guess which keep the problems associated with non-linearity 
under control. 

An effective way to construct leading approximations is to build Z in terms of spectral 
functions (propagators). A good combination can be guessed by writing down the per¬ 
turbation theory or an effective theory for the action being studied. However, the only 
information we use from these structures is the general spectral forms. The masses and 
couplings are initially arbitrary and are set by Galerkin. Another effective and rapidly 
convergent guess can be made using polynomials weighted with propagators. There are a 
huge number of reasonable choices. In practice, we have found remarkable stability (per 
computational time) between very different types of iterations schemes when used on simple 
problems. 

As a special case of the above non-linear form, we have been examining a spectral 
iteration scheme that we would think has much of the structure of an actual solution. The 
initial ansatz is taken to be roughly of the form 

\uZ = J dwdx-J{w)G 2 {w,x)J{x) 

+ J dwdxdydz-J{w)J{x)G2iw,x,y,z)J{y)J{z) (10) 

with 


and 


G*2ip) 
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( 11 ) 


Gl{p,q,r) = Gl{p)G^{q)j—^-^-^Gl{r)G^{p + q-r)^ (12) 

We are generating results based on this starting form in 4 dimensions on fairly big 
systems. We are finding that this method consumes small amounts of GPU (relatively) but 
at the cost of large amounts of memory. 


6 Boundary Conditions 

The discussion up to this point has avoided a fundamental issue. One of the most compli¬ 
cated and interesting problems associated with this type of numerical solution lies in the 
analysis of the appropriate boundary conditions imposed on the source differential equa¬ 
tions. If we look at the special case of zero space time dimensions (ultralocal model) there is 
only one differential equation. It is a third order differential equation so its solution involves 
the specihcation of three a priori arbitrary constants |^]. In particular, the normalization 
of Z, and the one and two field Green’s functions are unspecified. The situation for ar¬ 
bitrary space-time dimension is more complex with a correspondingly greater number of 
unspecified constants. This situation is not entirely unfamiliar. In this discussion we have 
excluded terms involving odd powers of the sources. This constitutes a boundary condition 




on the first derivative of Z. We know that, if we do not exclude odd terms, we can under 
the correct conditions write solutions that have spontaneous symmetry breaking. Dealing 
with the boundary condition on the second derivative terms is a much more complicated 
and a very interesting problem that leads to the enhanced understanding of phase structure 
mentioned earlier. In particular we find the possibility of solutions to quantum field theory 
which are extremely singular for small bare couplings. We analyze this elsewhere iSi- 
I have avoided direct contact with this problem in the example detailed here by truncating 
the Taylor series solution vn Z. It is easy to see that truncation forces the approximation 
to agree with the smooth g = 0 limit of the theory that corresponds to perturbation the¬ 
ory. Therefore, truncation imposes an implicit boundary condition on the second derivative 
terms. 


7 Conclusions 

Our previous published work in developing the Source Galerkin method has been focussed 
primarily on scalar in various spacetime dimensions and lower dimensional lattice four 
Fermi theories i i 0- Now that we have an understanding of these systems we have 
begun to analyze gauge and fermionic theories in four dimensions. We have developed more 
sophisticated expansion that the direct series methods and have a good understanding of 
how to use Galerkin methods to set parameters in effective held theories. We believe that 
Galerkin techniques can be effectively used on the continuum and have verihed this in two 
dimensional models. We have developed what appears to be a correct working procedure 
to deal with continuum four dimensional problems. 

As we develop our techniques we anticipate to be able to answer questions about more 
realistic theories. At this stage, I am very optimistic that the Source Galerkin model will 
be important in making a contribution to our understanding of the solutions of quantum 
held theory. 
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